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Abstract. Motor proteins are active enzyme molecules that play a crucial role in 
many biological processes. They transform the chemical energy into the mechanical 
work and move unidircctionally along rigid cytoskclcton filaments. Single-molecule 
experiments suggest that motor proteins, consisting of two motor domains, move in a 
hand-over-hand mechanism when each subunit changes between trailing and leading 
positions in alternating steps, and these subunits do not interact with each other. 
However, recent experiments on heterodimeric kinesins suggest that the motion of 
motor domains is not independent, but rather strongly coupled and coordinated, 
although the mechanism of these interactions are not known. We propose a simple 
discrete stochastic model to describe the dynamics of homodimcric and heterodimeric 
two-headed motor proteins. It is argued that interactions between motor domains 
modify free energy landscapes of each motor subunit, and motor proteins still move via 
the hand-over-hand mechanism but with different transitions rates. Our calculations 
of biophysical properties agree with experimental observations. Several ways to test 
the theoretical model are proposed. 

E-mail: tolya@rice.edu 



Interaction Between Motor Domains in Kinesins 



2 



1. Introduction 

Several classes of enzyme molecules that convert chemical energy into mechanical 
motion are called motor proteins, or molecular motors. In recent years, these 
proteins have attracted a significant attention because of their importance for multiple 
biological processes [U EJ [3j HI [5] . Motor proteins, such as kinesins, myosins, dyneins, 
polymerases and helicases, move in a linear fashion along the rigid biopolymers (actin 
filaments, microtubules, DNA and RNA molecules). Typically, a fuel for the motion 
of these nanomotors comes from the hydrolysis of adenosine triphosphate (ATP) or 
related compounds. Although some progress in understanding the mechanisms of 
motor proteins has been achieved [31 HI [5], there are still many unresolved issues. 
One of the most important fundamental questions concerning motor proteins is how 
different domains of these enzymes coordinate and regulate their complex dynamics 
and functions. The goal of this paper is to address some aspects of this issue from the 
theoretical point of view. 

The enzymatic activity of motor proteins takes place in the so-called motor subunits 
that contain ATP-binding sites. Motor proteins typically have several such domains. 
The functioning of molecular motors strongly depends on the relative position and 
dynamics of these subunits [31 [5]. Two possible mechanisms for two-headed motor 
proteins have been proposed: an inchworm motion and a hand-over-hand mechanism 
[3j [6] . In the inchworm mechanism one motor domain is always in the leading position, 
while another one always trails. However, in the hand-over-hand mechanism the motor 
domains alternate their leading and trailing positions as the motor protein molecule 
proceeds along the filament track. Single-molecule experiments that utilized fluorescent 
imaging with one-nanometer accuracy (FIONA) have shown that individual double- 
headed kinesins, myosins V and VI, and cytoplasmic dyneins step in the hand-over- 
hand fashion [TJ El El [101 E] • Thus, this mechanism explains stepping dynamics of the 
majority of motor protein species. 

In the current version of the hand-over-hand mechanism it is assumed that two 
heads move independently from each other, i.e., when the trailing motor subunit moves 
its dynamics is not affected by the presence of another motor subunit. Then a mean 
dwell time to advance one step forward for a heterodimeric motor protein with two 
heads labeled as A and B is given by 

ta/b = ^{t~a/a + tb/b), (1) 

where ta/a an d tb/b are mean dwell times for homodimeric A/ A and B/B motor 
proteins, respectively. The corresponding relation for the velocity can be written as 

ZVa/aVb/b , . 

V A/B = 77 —T7 • (2) 

VA/A + Vb/B 



However recent single-molecule investigations of dynamics of kinesins [T2] do not support 
these relations, and, consequently, the independence of two motor domains during the 
motion is put in a doubt. In these experiments force- velocity curves and enzymatic 



Interaction Between Motor Domains in Kinesins 



3 



activities have been measured for different homodimeric and heterodimeric kinesins. 
Surprisingly, it was shown that the velocity of the heterodimeric kinesin with a mutation 
in one of the motor heads is not given by Eq. (2). It was suggested that the two heads 
strongly influence the dynamic and enzymatic properties of each other, although the 
mechanism was not specified. In this paper we present a simple discrete stochastic model 
that might explain several aspects of the complex dynamics of heterodimeric kinesins. 
Our main idea is that motor domains interact with each other and significantly modify 
the overall dynamics. 

Theoretical investigations of molecular motors follow several approaches that 
include continuum ratchets [13], [TJ], discrete stochastic models [5] and computer 
simulations [151 HE]- I n this work we utilize a discrete stochastic approach because 
it conveniently provides explicit expressions for dynamic properties, and it is able 
to describe successfully different aspects and trends of motor proteins transport 

2. Theoretical Model 

In our theoretical model it is assumed that a kinesin protein molecule steps from a 
given binding site to the next one at a distance d = 8.2 nm along the microtubule 
through N intermediate biochemical states. Different kinetic schemes for the motion of 
homodimeric and heterodimeric kinesin molecules are shown in Fig. 1. Similar approach 
has been used successfully to describe dynamic properties of kinesins and myosins V 

EH eh mi- 
Let us label homodimeric motor proteins of wild type as W/W, homodimeric motor 
proteins with mutations in both heads as M/M, and heterodimeric proteins with only 
one mutated motor domain as W/M. Although the motion of motor proteins includes 
many complex biochemical and biophysical processes, in the simplest approximation it 
is reasonable to use only two-state (N = 2) discrete stochastic models to describe the 
dynamics of both homodimeric kinesins: see Fig. 1. The motor protein molecule can 
jump forward from a state to a state 1 with the rate uq, and this transition corresponds 
to ATP binding, yielding 

«o = A;o[ATP], (3) 

where ko is a rate constant. The reverse transition is given by the rate w±. The forward 
and backward transitions between the state 1 and the state on the next binding 
site, with the rates u\ and wo respectively, describe several biochemical processes, such 
as ATP hydrolysis and release of hydrolysis products, but combine them in one step. 
We assume that mutations do not affect strongly ATP binding process, but only the 
enzymatic functions are changed, which leads to a different pair of transition rates (u[ 
and w'q) for the mutated molecule (see Fig. 1). Experiments show that some mutations 
decrease the binding affinity of motor proteins to microtubules, thus leading to decrease 
in enzymatic activity |12j . The transition rates for W/W and M/M motor proteins are 
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Figure 1. General schematic view of discrete stochastic models for kinesins. W/W 
labels the homodimcric motor protein with both wild-type motor heads, M/M 
corresponds to the homodimcric kinesins with mutations in both motor domains, while 
W/M describes the heterodimeric motor proteins with mutation in only of the heads. 
Parameters e, £\ and £2 describe free energy changes due to mutations relative to the 
wild-type motor proteins. 



related via the detailed balance condition, 

Ui 



u 



Wn 



exp(-e/k B T). 



(4) 



^0 w o 

The parameter e describes the effect of mutation on the enzymatic properties of 
the motor protein, i.e., how the ATP hydrolysis for the mutated homodimer is 
thermodynamically less favorable in comparison with that of wild-type homodimer. 
Microscopic origin of this parameter is the result of complex interactions between two 
motor domains and between motor heads and the microtubule track. 

The situation is more complex for heterodimeric W/M kinesins, as illustrated in 
Fig. 1. Two motor domains and interactions between them are different from the cases 
of homodimer kinesins. As a result we have two different sets of rates to model ATP 
hydrolysis by the wild head (u^ and w^) and by the mutated head (uf' and w^). 
These rates are also related via detailed balance conditions, 



(i) 



w 



(1) 



Wq 



exp(-e 1 /k B T), 
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(2) 



It is important to note that generally, E\ ^ Ei 7^ e, because of different interactions 
between the motor domains. Then, assuming that the hand-over-hand mechanism still 



(2) 



Wq 



exp(-e 2 /k B T). 



(5) 
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is a valid description for the stepping of individual molecules, the dynamics of W/M 
kinesin molecules can be described by iV = 4-state model with step size equal to 2d=lQA 
nm. 

The explicit expressions for the rates can be obtained from Eqs. fl4]) and ([5]): 



W17 



u 



«i7i 



(2) 



Ml7 2 



02 . 



and 



w, 



o 



w 7 



l-a 



(1) 



w 7i 



,1— Cfl 



(2) 



where we defined 

7 = exp 



, 7i = exp 



72 = exp 



f-V 



(6) 
(7) 
(8) 



\kT) 

Parameters a, «i and a 2 are energy-distribution factors that describe how the free 
energies change due to the mutation affects corresponding to the forward and the 
backward transitions. For simplicity, we assume that a = a\ = a 2 , although more 
general situations can be easily analyzed. 

The advantage of using discrete stochastic models is the fact that all stationary- 
state dynamic properties of motor proteins, such as mean velocities, mean dispersions 
and stall forces, can be obtained exactly for any number of intermediate states N in 
terms of the forward (uj) and backward (wj) transition rates (5j [19], [21] . Specifically, 
exact expressions for the mean velocity can be presented in the following form, 



K N 



N-l 

i-TT^ 



j=0 J 



(9) 



where d is the step size (equal to 8.2 nm for N = 2 model and 16.4 nm for N = 4 
model), and the auxiliary functions Rn are given by 



R 



N 



T 3, T i 



1 

U ; 



N-l 

3=0 3 

Specifically, for the velocity of homodimeric motor proteins we obtain 



N-l j+k 

En? 

k=l i=j+l 



U + Ui + W + Wi 



M Ml7~ 



w wa 



l-a 



V(W/W) = d 
for W/W kinesins, and 

V(M/M) = d 

L u o + ^i7" a + ^o7 
for M/M kinesins. For heterodimeric motor proteins Eq. (Q yields 

"1^(7172)"" - u'o U7 i(7i72) 1 " Q " 



V(W/M) = 2d 



A 



with the parameter A given by 

A = (li a + l2 a )u ui(u + Wl ) + (71 1 " + 72 _a )^o^i(wo + wi) + 

(7i72)~°(7i + 72)^1^0(^0 + wi) + 2(7i72)~ q (moM? + WqWi7i72)- 



(10) 



(12) 



(13) 



(14) 
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The general expression for dispersion in sequential discrete stochastic models can 
be written in the following form [5], [191 12T] - 

D = (d/N) {[VS N + dU N ]/(R N ) 2 — (N + 2)V/2} , (15) 

where 

N-l N-l N-l 



j=0 k=0 j=0 J 



N-l j-k 

i+£ n 



Uj 

k=l i=j-l 1 



The explicit equations for dispersions of W/W, M/M and W/M kinesins can be obtained 
in the way similar to the velocities, however, these expressions are quite bulky and they 
will not be presented here. 

When the motor protein is a subject of external loads, the resisting force that 
completely stops the molecule is called a stall force, F$. For general iV-state sequential 
discrete-stochastic models the stall force can be written as [5j [17] 

For homodimeric kinesins our model predicts the following stall forces, 
F S (W/W) = —— In , Fs{M/M) = —— In 



d wqWi d wqWi^ 

Comparing these equations, we obtain 

F S {M/M) = F S (W/W) - e/d. (19) 

For heterodimeric kinesins the stall force is given by 



k B T ulu\ 
2d wlw\ 7i7 2 : 



which leads to 

F S {W/M) = F S (W/W) - (e x + e 2 )/2d. (21) 

Eqs. (I19p and ([2~TT) provide a simple physical interpretation and a method of estimating 
parameters e, Z\ and e 2 . 

External force F also strongly modifies transitions rates [3 [T7] 

""fcr )' (22) 

6rFd\ 

Wj (F)=^(0)expl+-^ T j, (23) 

where 6 • are load-distribution factors that describe how the external load changes the 
energy activation barriers for the forward and backward biochemical transitions from 
the state j. Load-distribution factors are related via 

N-l 

£W + *7) = 1 - (24) 

3=0 



.(16) 
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3. Results and Discussion 

In the experimental work of Kaseda et al. [12] the coordination of two heads for 
different homodimeric and heterodimeric kinesin molecules has been investigated using 
microtubule-gliding assays and optical trapping spectroscopy. Different homodimeric 
and heterodimeric motor proteins were prepared by mutations in the motor domains 
that affect microtubule-binding region. It was found that dynamic properties of 
heterodimeric proteins with one mutated head could not be described by independent 
hand-over-hand stepping mechanism. 

To analyze experimental data we consider kinesins with only one type of mutation, 
although our method can be easily applied to different molecular motor species. 
Specifically, it was shown [12] that at [ATP] = 1 mM a wild- type homodimer travels 
with the stationary velocity V(W/W) = 679 ± 59 nm/s, and it produces the maximum 
stall force of F S (W/W) = 6.3 ± 0.9 pN. When the mutation labeled L12 affects both 
motor heads the resulting homodimeric M/M kinesin does not attach to microtubules, 
indicating zero velocity and stall force. However, surprisingly, heterodimer with the same 
mutation L12 in one of the motor heads can move with the velocity V(W/M) = 101 ±25 
nm/s, while exerting the maximal stall force of Fg(W/M) = 0.8 ± 0.2 pN. Then, from 
Eqs. (HI and (J2U) we obtain 

e = 12.6 ± 1.8 kgT, e Y + e 2 = 11 ± 4.4 k B T. (25) 

The important result is that E\ + e 2 < 2e, which indicates that biochemical properties of 
mutated and wild-type motor heads in the W/M kinesin differ from the corresponding 
properties in the W/W and M/M motor proteins, supporting our idea of modifying 
free-energy landscapes for motor proteins via interaction between motor subunits. 

After systematically exploring parameter space, we found that all experimental 
data for kinesins with L12 mutation can be well described by the following parameters, 

k = 1.2^M _1 s -1 , Mi = 90s" 1 , w 1 = 10s~\ w = 0.05s" 1 , 

a = 0.14, e = 12.3 kgT, e x = e 2 = 10.6 kgT. (26) 

Note that obtained transition rates are similar to parameters previously utilized to 
describe the dynamics of single kinesins, and they are consistent with chemical kinetic 
experimental results [21] . In addition, in our calculations we also used load-distribution 
factors reported earlier [21] . 

0+ = 0.135, 07/ = 0.080, 0+ = 0.035, 9 = 0.750. (27) 

The results of theoretical calculations for the velocities of homodimeric and 
heterodimeric kinesins at different conditions are presented in Fig. 2. The effect of 
external load on the motor protein dynamics is shown in Fig. 2a. The force slows down 
the motion of all motor proteins, as expected. It is found that the stall forces Fg for 
W/W and W/M kinesins are equal to 6.2 pN and 0.8 pN respectively, which are in 
excellent agreement with experimentally measured values [12] . The force- velocity curve 
for heterodimer is essentially linear, while for homodimeric W/W kinesins it deviates 
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from the linear dependence. Again, this theoretical result agrees well with experimental 
observations (compare with Fig. 4b from Ref. [12] ). The dependence of the velocity on 
the concentration of ATP (see Fig. 2b) shows typical Michaelis-Menten character, i.e., 
a linear behavior at small concentrations of ATP and a saturation at large [ATP]. One 
of the important theoretical predictions of our model is that the heterodimeric kinesin 
with L12 mutation in one of the motor heads will not function (V = 0) even without 
external resisting force (F = 0) for ATP concentrations smaller than 0.186 mM. This 
theoretical prediction can be easily tested in experiments. 

In addition to the mean velocity, an important characteristic of the motor protein 
dynamics is a diffusion constant or dispersion [5]. In our theoretical framework it can 
be calculated exactly, and the results of these computations are given in Fig. 3. The 
effect of external forces on dispersions for homodimeric and heterodimeric kinesins is 
shown Fig. 3a. External forces have qualitatively different effects on the fluctuations 
of wild- type and mutated motor proteins. Dispersion for homodimeric W/W kinesins 
is a decreasing function of the external load, while for homodimeric M/M proteins it 
increases. It is interesting to note that the behavior of the heterodimeric W/M is non- 
monotonic, depending on the strength of the external forces. At low external loads 
(F < 1 pN) the dispersion of W/M kinesins decreases, as is for W/W kinesins, however 
for larger external forces it starts to increase, similarly to M/M kinesins. Dispersion as 
a function of [ATP] is plotted in Fig. 3b. Again, the behavior of heterodimeric W/M 
proteins is intermediate between wild-type and fully mutated homodimeric kinesins. 

Our theoretical calculations of dynamic properties and analysis of the experimental 
observations for homodimeric and heterodimeric kinesins suggest the following 
mechanistic (however strongly simplified) picture of the motor protein dynamics. 
Mutations change the free energy landscapes for the enzymatic activity and the 
mechanical processivity of molecular motors. However, motor domains in W/W, M/M 
and W/M proteins interact with each other differently, leading to different free energy 
surfaces. As a result, two heads in the heterodimeric molecule become very similar 
in the biochemical properties, but different from the corresponding motor domains in 
the homodimers. The heterodimeric motor protein still moves along the microtubules 
in the hand-over-hand mechanism, although with the transition rates modified by 
interaction between motor heads. We suggest that this theoretical picture can be tested 
in experiments by, for example, labeling differently both motor domains to obtain a 
necessary dynamic information. 

4. Conclusions 

We developed a simple theoretical description of the dynamics of motor proteins based 
on the discrete sequential stochastic models. This approach allows us to resolve the 
contradiction between experimental observations on homodimeric and heterodimeric 
kinesins and the widely accepted hand-over-hand stepping mechanism for two-headed 
molecular motors. It is argued that the interaction between motor domains can modify 
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Figure 2. a) Force- velocity curves for W/W homodimeric kinesins (black line) 
and for W/M heterodimeric kinesins (red line) at [ATP] = 1 mM. b) Velocities of 
W/W homodimeric kinesins (black line) and W/M heterodimeric kinesins (red line) 
as a function of [ATP] at the constant external force F = 0.5 pN. The velocity of 
homodimeric M/M kinesins is zero at all conditions. 
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Figure 3. a) Dispersion as a function of the external force for W/W homodimeric 
kinesins (black line), for homodimeric M/M kinesins (green line), and for W/M 
heterodimeric kinesins (red line) at [ATP]=1 mM. b) Dispersions of W/W homodimeric 
kinesins (black line), homodimeric M/M kinesins (green line), and W/M heterodimeric 
kinesins (red line) as a function of [ATP] at the constant external force F = 0.5 pN. 
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free energy landscapes for the motor protein motion, and the transitions rates change 
depending on the nature of these domains. Explicit calculations of dynamics properties, 
such as velocities, dispersions, and stall forces are presented for homodimeric and 
heterodimeric kinesins with L12 mutations. Theoretical predictions agree well with 
available experimental data. Several suggestions on how to test theoretical predictions 
are discussed. 
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